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Abstract 

We present a study of the normal state optical conductivity in the cuprate 
superconductors using the nearly antiferromagnetic Fermi liquid (NAFL) de- 
scription of the magnetic interaction between their planar quasiparticles. We 
find that the highly anisotropic scattering rate in different regions of the 
Brillouin zone, both as a function of frequency and temperature, a bench- 
mark of NAFL theory, leads to an average relaxation rate of the Marginal 
Fermi Liquid form for overdoped and optimally doped systems, as well as for 
underdoped systems at high temperatures. We carry out numerical calcula- 
tions of the optical conductivity for several compounds for which the input 
spin fluctuation parameters are known. Our results, which are in agreement 
with experiment on both overdoped and optimally doped systems, show that 
NAFL theory explains the anomalous optical behavior found in these cuprate 
superconductors. 
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I. INTRODUCTION 



Optical conductivity measurements show that high temperature superconductors exhibit 
a number of anomalies when compared to the usual Drude-like Fermi liquid (FL) behavior 
found in conventional metals. Quite generally, experiment shows that the effective trans- 
port scattering rate, 1/t(u), exhibits linear in uj behavior over a wide frequency range, 
while the measurements of Puchkov et al! demonstrate that significant deviations from this 
standard behavior occur, especially in underdoped materials at low temperatures, in the so 
called pseudogap regime, where 1/t{uj) is strongly suppressed. These anomalies are often 
attributed to a possible non-Fermi liquid ground state in these materials.! 

The linear frequency dependence of 1/r at relatively high u can be regarded as the 
optical analogue of the anomalous resistivity, p oc T, found in the superconducting cuprates, 
and its presence has been considered as major support for the marginal-Fermi liquid (MFL) 
approach.! In the MFL approach, the self energy is given by: 



where x = max(|u;|, T), u c is a cutoff frequency and A a coupling constant, an expression 



spectral weight and lifetime are highly anisotropic as one goes around the Fermi surface 
(FS), rather than being momentum independent as Eq. (JJ) requires. 

Highly anisotropic quasiparticle behavior finds a natural explanation in the nearly an- 
tiferromagnetic Fermi liquid (NAFL) model of cuprates.! in which the effective interaction 
between planar quasiparticles is assumed to be proportional to the dynamic spin suscepti- 
bility, x(g,co>), measured in NMR and INS experiments. Because u;) is strongly peaked 
in the vicinity of the antiferromagnetic wavevector Q = (tt,tt), quasiparticles on the FS 
located in the vicinity of the hot spots,! portions of the FS which can be connected by 
wavevectors ~ Q, interact strongly, while those located elsewhere (in cold regions) interact 
comparatively weakly. 




(1) 
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NAFL theory, with spin-fluctuation parameters and a quasiparticle spectrum taken from 
experiment, yields results consistent with low frequency NMR and INS experiments, as well 
as ARPES measurements of the Fermi surface properties.0! We have recently shown that 
it also provides a consistent qualitative, and in many cases quantitative, explanation of the 
measured changes with doping and temperature of both the longitudinal and Hall conduc- 
tivities of the cuprates, and shown how one can extract the quite distinct lifetimes of "hot" 
and "cold" quasiparticles as a function of temperature directly from these measurement s.@ 
In this paper we consider NAFL theory at infrared frequencies and show that it explains as 
well the measured anomalous optical behavior. 

Our paper is organized as follows. Following a review of the physical properties of NAFLs, 



we discuss in Sec. Ill both analytic and numerical calculations of the highly anisotropic quasi- 



particle lifetime and the quasiparticle relaxation rate at finite frequency u and temperature 



T. In Sec. [Ty| we calculate the optical conductivity and discuss its limiting behavior. In 
Sec. [V] we use realistic band and spin-fluctuation parameters to compare our results with 
experiment. Our conclusions are presented in Sec. |VJ. 

II. THE PHYSICAL PROPERTIES OF NAFLS 

In the NAFL description of the normal state properties of the superconducting cuprates, 
it is the magnetic interaction between planar quasiparticles which is responsible for the 
anomalous spin and charge behavior. The magnetic properties of the system are specified 
by the dynamical spin-spin response function of fermionic origin, x(q, u), which near a peak 
at a wavevector Qj in the vicinity of Q, is assumed to take the mean-field form:0 

x(q,u) = —-, — r^cz — ~i — ( 2 ) 

1 + (q - Qi) 2 <r - ioj/uJ s f 

Hence \Q — a £ 2 ^ Xo is t ne magnitude of the static spin susceptibility at Qj, £ is the anti- 
ferromagnetic correlation length, u 8 f specifies the low frequency relaxational mode, brought 
about by the near approach to antiferromagnetism and a is a temperature independent 
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scale factor. We use a system of units in which the lattice spacing a = 1. Although ex- 
periment suggests that the susceptibility (0) quite generally possesses four peaks at the 
incommensurate wavevectors Qj ~ Q, since the introduction of the incommensuration does 
not introduce qualitatively new physics we shall assume for the most part that the spin 
fluctuation spectrum possesses only a single peak at Q. 

The quasiparticle spectrum is assumed to take a tight-binding form, 

6k = — 2t(cos k x + cos k y ) — At' cos k x cos k y — 2t" (cos(2k x ) + cos(2k y )) (3) 

where t, t' and t" are the appropriate hopping integrals, while the effective magnetic inter- 
action between the planar quasiparticles is specified by 

V eff (q,u) = g 2 x(q,u). (4) 

For a given system the parameters Xq,h 6; ^sf which determine x(q, u) are taken from fits 
to NMR and INS (inelastic neutron scattering) experiments, while the effective coupling 
constant, g, is assumed to be momentum independent for wavevectors near Qj. As dis- 
cussed by Chubukov et al,0 the effective interaction, Eq. (f|), can be, in principle, derived 
microscopically starting with, e.g., a one-band Hubbard model. 

The effective interaction V e ff, Eq. (f|), has the obvious property that for sufficiently 
large correlation lengths it is highly peaked for momentum transfers in the vicinity of the 
antiferromagnetic wavevector Q. The consequences of this peaking are hard to overestimate: 
if the FS of the system of fermions, defined by the quasiparticle dispersion (|3|), is such that 
it intersects the magnetic Brillouin zone (see Fig. [l|), then quasiparticles in the vicinity of 
these intersection points on the FS, the hot spots,! are much more strongly scattered by the 
spin-fluctuations than those which are on other parts (cold regions) of the FS. This leads to 
strongly anisotropic quasiparticle behavior, since the resulting temperature and frequency 
variation of the quasiparticle scattering rates at and far away from hot spots is in general 
very different. We shall return to this point in the following Section. 

Barzykin and Pines (BP) have used Eq. (Q) to analyse NMR results in the cuprates.0 
They find that for underdoped systems the low frequency magnetic behavior possesses three 
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distinct normal state phases, and that the characteristic frequencies and lengths which enter 
Eq. (Q) are connected by a dynamic scaling relationship, u s f oc £~ z , in two of the phases. 
Thus, above the temperature, T cr , at which the temperature dependent uniform suscep- 
tibility, xo(T), possesses a maximum, the system exhibits non- universal mean field (MF) 
behavior with dynamical exponent z = 2; u s f ~ l/£ 2 , and the product, XQ^sf ~ u sf£, 2 is 
independent of temperature. From a detailed analysis of the NMR experiments, Barzykin 
and Pines conclude that the crossover temperature, T cr , is determined by the strength of the 
AF correlations, with £(T cr ) « 2. Below T cr , down to a second crossover temperature, T*, 
underdoped systems exhibit non- universal, z — 1, pseudoscaling (PS) behavior. Thus for 
T* < T < T cr , it is u) s fl; which is independent of temperature. NMR experiments show that 
above T cr , in the MF regime, u s f and l/£ 2 scale linearly with T, i.e., u s f = A + BT, while 
in the PS regime, between and T cr , it is u s f and which scale linearly with T, albeit 
with a possibly somewhat different slope and intercept of uj s f than that found above T cr . 
This behavior has now been verified experimentally for the YBa2Cu307_ :E , La2_ z Sr x .Cu04 
and YBa2Cu40g systems,0 as well as the Hg and 2212 BSCCO compounds.0 Below T*, 
in the pseudogap (PG) regime, £ becomes independent of temperature while u> s f, after ex- 
hibiting a minimum near T*, rapidly increases (roughly as 1/T) as T decreases toward T c . 
These changes in the magnetic fluctuation spectrum are accompanied by (indeed result from) 
changes in the quasiparticle spectrum.0 Thus, the PS regime is characterized by a strong 
temperature variation of the quasiparticle spectrum, resulting in a FS evolution0 which has 
non-trivial consequences for the transport properties of underdoped systems.! The crossovers 
seen in the low frequency magnetic behavior thus possess, to a considerable extent, their 
charge counterparts in transport experiments. 

From a magnetic perspective, the so called optimally doped systems (e.g., YBa 2 Cu 3 6 .93 
and Lai. 85 Sro.i5Cu0 4 ) are a special case of the underdoped systems, in which T* is com- 
paratively close to T c . Overdoped cuprates are defined as those for which T cr < T c . For 
these systems, then, the antiferromagnetic correlations are comparatively weak, with £ < 2, 
Xo{T) is at most weakly temperature dependent, while u s f oc £ -2 follows the linear in T 
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behavior found in the underdoped systems above T cr . Examples of over doped systems are 
the T c ~ 40K Tl 2212 system and La 2 _ x Sr 1 .Cu0 4 for x > 0.24. 

As Monthoux and Pines,0 and Chubukov et al0 have emphasized, the physical origin 
of the highly anomalous behavior displayed by the nearly antiferromagnetic Fermi liquids 
resides in the non-linear feedback of changes in the hot quasiparticle spectrum on the strong 
quasiparticle interaction, ~ x(q, u), which determines that hot quasiparticle behavior. For 
systems above T cr , or in overdoped systems, that feedback is negative; strong coupling ef- 
fects cause Xo{T) to increase weakly as the temperature decreases, while an RPA or mean 
field description description suffices to determine the relationship between u 8 f and £. At a 
critical value of the strength of the AF correlations, that feedback becomes positive; a weak 
pseudogap develops in the hot quasiparticle spectrum, bringing about a linear decrease in 
Xo(T) as T decreases below T cr , while the quasiparticle damping of spin excitations becomes 
T dependent in such a way that u s f displays z = 1 pseudoscaling behavior, with the strength 
of the AF correlations, £, growing as (a + bT)^ 1 . The second crossover temperature, T*, 
marks the transition to "strong" pseudogap behavior: the AF correlations become frozen; 
u s f increases rapidly as T decreases below T*, while Xo(T) falls off rapidly between T* and 
T c . It is the interplay between these changes in the quasiparticle spectrum (seen directly in 
ARPES measurements) and the spin fluctuation spectrum, and effective quasiparticle inter- 
action, seen directly in NMR and INS experiments, which is responsible for the anomalous 
quasiparticle lifetimes we now consider. 



III. QUASIPARTICLE LIFETIME AT FINITE FREQUENCY 

We begin by estimating the effective scattering rate, for quasiparticles near the FS, as a 
function of frequency and temperature using simple perturbation theory. The quasiparticle 
lifetime is determined by the imaginary part of the single particle self-energy, which, for the 
effective interaction in the second order Born approximation, reads: 
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ImZ(k, uj) = Y,9 2 Imx(k -U,u- e') [n (e' - u) + / (e')] (5) 

it' 

where e' = e^, and n (e), and /o( e ) are the Bose and Fermi distribution functions respec- 
tively. The transport relaxation rate is usually given by: 
1 r d 2 k' 



J -^2g 2 Im X (k - k', d - u) [n (e' - u>) + f (e')](l - lk , k .) (6) 



T k (u) J (27r) 

where is a temperature independent vertex function which removes forward scattering 
processes from the calculated transport relaxation rate. As we shall see below, for overdoped, 
optimally doped, and even underdoped systems above T*, the dominant contribution to 
l/rfc comes from large momentum transfer scattering processes and hence the use of 7^/ 
for forward scattering corrections plays only a minor role. Therefore we assume 7 = 0, and 
return later to the role of the vertex corrections at low temperatures in underdoped systems. 
We make the usual change of variables in the integral in Eq. (0): 

where k' is the component of momentum k' parallel to the equipotential lines (e' = const), 
and carry out the integration over e' for k near the FS and e < t, where t is the hopping 
matrix element. On making the substitution x = (V — u))/T in Eq. (^|), it becomes: 

where 



2ag 2 u sf e [ fV(k,k',e) (8) 

J FS 47T Z 



I(k,k , ,e) = J 



dx 



1 1 

+ 



exp(x) — 1 exp(x + e) + 1 



x 



n 2 + x 2 



(9) 



and 



e = uj/T fi = ^ / [l + e 2 (k-k , -Q) 2 ]/T = ^/T. (10) 

Note that the definition ([TDD effectively removes the temperature from the calculation. 

Obviously, the integral / is hard to solve analytically. However we can obtain its value 
in appropriate limits, and then map its behavior with a universal function which possesses 
the correct limiting values. 



Low Frequency Limit, e <C 1: In this limit one can expand l/[exp(x + e) + 1] in Eq. @. 
The first non-vanishing correction to I(e = 0) is l2£ 2 , where 

I y°° x d 2 ( 1 \ 

2 2 Loo x 2 + ft 2 &r 2 ^exp(x) + I) ' 1 J 

When f2 -C 1 one finds I2 — 7C(3)/27r 2 — 7rfi/8 to lowest order in Q, while in the limit Q ^> 1 
one finds /2 ~ I/O 2 to lowest order in I/O. 

Low Temperature Limit, e > 1, fi > 1: In this limit, to leading order in 1/e and 1/Q 
Eq. (|D becomes: 

1 / e 2 \ 7T 2 Q 2 - e 2 n 2 

^) = 5 lo 4 1 + s?j + -6-(sFTP? + 3S5 (12) 

The first term in (|l|) is the zero temperature result which correctly incorporates the above 
small e behavior. The remaining two terms are of order T 2 , as in usual FLs. If the temper- 
ature is increased while e 3> 1 is maintained, then the result depends on Q in a nontrivial 
way: a quick check reveals that in the Q <C 1 limit / « tt/Q, as in the zero frequency case. 

If one seeks a universal function for I, it is not difficult to see that an appropriate form, 
for all relevant values of Q and e, is: 

'M = 5 log ( 1 + S + 2H(ST^ < 13 > 

This function displays the correct limiting behavior for both low and high frequencies and 
temperatures. As a function of Q it has a crossover from 1/Q to l/f2 2 behavior, reflecting 
the fact that as a function of temperature / displays a crossover from {-kT) 2 / uo\ k , at low 
temperatures to nT/ukk' at high temperatures (see Eq. ([TOj)) . 

The result (|TJ|) leads to the correct behavior of the scattering rate in the case of typical 
metallic conductors, for which the relevant energy scale u^y is of order Fermi energy, Ef, 
Uk^y is weakly momentum dependent and Q ~ Ef/T S> 1 for all k and k' . On using Eqs. 



(O) and (p~3|) , one trivially recovers the FL result: 



1 uo 2 + (vrT) 2 



oc (14) 



'f 



It is important to realize that the form of the integral / in Eq. (|T3| ) is such that one can 
treat the dependence on frequency and on temperature almost independently. When e and 
Q are of comparable size, the second term in Eq. fll^D is small compared to the third, and 
finds a result, which, although differing somewhat from Eq. fllcf), 



(15) 



turns out to be correct for a wide variety of experimentally relevant parameters. For Eq. 
(|15|) the crossover from n/fl behavior at Q <C 1 to tc 2 /Q 2 behavior at Q ^> 1 occurs at 
considerably lower values of Q than is the case with Eq. (|T"3|) . 



For the case at hand, for which uj^y < kT, we follow Refs. £)] and [15] and carry out the 
integration over k! in Eq. (B) by making the approximation 



UJ kk , W LU s f 



(16) 



where Ak' and Ak measure the displacements of k and k! from the corresponding adjoint 
hot spots along the FS. Then, for doping levels such that one has a large FS, to a high 
degree of numerical accuracy, the integral over Ak! can be extended to infinity, and for a 
general point near the FS the scattering rate is given by: 

1 ag 2 iJUT s 



2%v f 



{*t + F > 



(17) 



where 



irT + u sf (l + ¥(Aky 



(18) 



and 



The behavior of F^, which determines the low frequency, uj independent limit of l/r/-( 



has been discussed in some detail in Ref. 0; we turn therefore to the behavior of F% in various 



limiting cases. 



Finite (but low) frequencies, u <C ou s f' As indicated by the e — > behavior of integral 
I{k, k', e) above, the low frequency correction to the (temperature induced) scattering rate is 
always proportional to e 2 oc uj 2 , as in ordinary FLs. However, unlike Landau Fermi liquids, 
from Eqs. ([19]) one finds that the coefficient of proportionality, is strongly momentum 
dependent: 

l/r k (T, uj) - l/r k (T, uj = 0)^ — ^/ U ^t> 2 2^ ^ 

For hot quasiparticles with Ak£ <C 1, the characteristic energy scale which determines this 
FL-like behavior is small, being ~ uj s f, while for cold quasiparticles, with A/c£ 3> 1, it 
is much larger ~ u s f^ 2 (Ak) 2 . Hence, in the cold regions 1/r is proportional to uj 2 up to 
rather high frequencies of order w s j-( 2 ~ 100 meV. This result has important consequences: 
although the range of frequencies where the FL-like scattering rate persists is very small in 
the hot regions, one can use FL theory to calculate the effective mass, m* = m(l — d'E'/doj), 
which may be large, but is always finite. For cold quasiparticles one finds a comparatively 
modest effective mass enhancement at all temperatures of interest. 

Intermediate frequencies, uj s f <C uj <C ^sf£, 2 '- It follows from Eq. ( [3~9|) that the hot 
quasiparticles possess a lifetime: 
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- « 7 1 m + ^~ ^u sf e (21) 

T k T k (UJ = 0) 2-KVf V 

while for cold quasiparticles, as long as uj <C uj s f^ 2 (Ak) 2 , Eq. (|20|) still holds. 

Intermediate frequencies, uj ~ uj s fC, 2 (Ak) 2 , Ak < 1: From Eq. fll9|) it is also relatively 
easy to obtain a simple expression for the cold quasiparticle lifetime at intermediate values 
of frequency and Ak, such that uj ~ uj s f^ 2 {Ak) 2 : 

1 \ « (22) 

Tfc Tfc (cj = 0) 2-KVfAk y ' 

where rj is a numerical coefficient of order 0.7. On combining Eqs. ( pO"D and (|22|) and making 
use of our previous results for l/r fe (c<j = 0), we see that, to a good degree of accuracy, we 
can write: 
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r k (uj,T) ~ Av f Ak \ttT + nT (k) + to + E (k) 



1 / (vrT) 2 



(23) 



where 



7rTo(fc) ~ E (k) ~ 2w s/ e 2 (Afc) 2 



(24) 



The quantities 7rT and _E are not identical, but are nevertheless of comparable magnitude. 
Hence we see, as might expected for any interacting FL, that the same energy scale Eo(k) 
(up to a multiplicative constant of order tt) determines the crossover from FL behavior at 
low uj and T to a non-Fermi liquid linear in u and T regime. The crossover frequency 
increases somewhat with increasing temperature if, as is the case in overdoped materials, 
uj s f no longer plays a negligible role in determining l/r^w). In this case E is large and 1/r 
is only weakly frequency dependent. 

As in the case of the d.c. resistivity, one expects the dominant scattering rates for the 
optical conductivity to be those of the relatively well defined quasiparticles in the cold regions 
of the BZ; the corresponding values of Ak lie within a range Ak ~ 1 — 2 for an overdoped 
material.© For an underdoped material Ak is much smaller, Ak < 0.5, since the FS is not 
only close to the magnetic BZ boundary, but there may be strong incommensuration effects. 
In addition, while in the overdoped materials w s /£ 2 = const as required by the z = 2 
scaling, in underdoped materials one encounters the z = 1 scaling relationship between the 
characteristic energy and momenta, where w s j( = const and £ ~ 1/T. 

A key feature of the NAFL model is that the relevant scale for the frequency and temper- 
ature variation of the relaxation rates, Eq = uj s f^ 2 {Ak)1 nax /^ is considerably smaller than 
the fermionic bandwidth. Hence at practically any frequency of interest (say, uj > lOmeV) 
and in the temperature regime where the resistivity is linear in T one finds from Eq. flTTp 
that the momentum average scattering rate as a function of frequency and temperature takes 
the MFL form: 




(25) 
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Thus, the NAFL approach leads to the major result of MFL theory^ without resorting to a 
model in which (see Eq. (|1|)) one encounters an infinite effective mass of the quasiparticles; 
indeed, as noted above, in the NAFL model the effective mass of the cold quasiparticles is 
always of the order of unity. 

Eq. ( |2"3"D implies a parallel between the T and uj dependence of the scattering rates: since 
the same characteristic frequencies appear in both cases^l the behavior of 1/r as a function 
of uj or T is approximately the same. We verify this in Figure |2| where we plot the zero 
temperature scattering rate as a function of uj at fixed large £ and uj s f — > for several 
values of Ak near the FS. The top panel shows quasiparticle scattering rates when the FS 
is large (Eq. fl2T|)), while the bottom panel displays our result for a small effective FS with 
a maximum distance from a singular point along the FS of Ak max = 0.3. In the latter case 
the integral over k' in Eq. (|8]) was done numerically. The corresponding scattering rates as a 
function of temperature at zero frequency are given in the insets of the panels. The parallel 
between the zero-frequency temperature and zero-temperature frequency dependencies of 
the scattering rate is quite obvious. Clearly, the rate is highly anisotropic as a function 
of momentum and decreases dramatically as one departs from the singular point, Ak = 0, 
where it is the largest and most anomalous, 1/r ~ y/u. On the other hand we notice that 
the influence of the FS size is felt mostly in the hot regions while cold quasiparticle lifetimes 
are comparatively insensitive to the size of the FS. 

We consider next the influence of changes in the spin fluctuations spectrum on the 
quasiparticle scattering rate. We begin with systems and temperatures which display z = 2 
scaling, where uJ s f(, 2 does not depend on temperature. As explained in Sec. ||, this scaling 
law is observed in NMR measurements on overdoped systems or underdoped systems at 
temperatures above T cr . In Figure 0a we show the frequency dependence of the hot quasi- 
particle scattering rate for several different temperatures. In order to make the connection 
with experiment, in this, and many subsequent figures we plot the inverse mean free path 
in the units of lattice spacing, by assuming a constant Fermi velocity, Vf = 0.25eV. As 
in previous work,! we have assumed realistic spin-fluctuation parameters^ appropriate for 
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YBa 2 Cu 3 7 , u sf = [6 + 0.06 T (K)} meV and oo sf ^ 2 = 60 meV. We see that the shape of 
the curve is somewhat less anomalous than that shown in Fig. |2|. Since the onset of the 
non-FL frequency dependence is at u ~ u s f for a hot spot, the observed behavior is a 
natural consequence of using a finite value of u s f and a correspondingly modest value of 
£. Since ui s f increases with temperature the scattering rate becomes even more FL-like at 
higher temperatures . Quite generally, outside of the pseudogap regime, the behavior of hot 
quasiparticles becomes less anomalous as T increases, since the system is further away from 
the AF instability, as measured directly by the size of l/£ and/or u s f. Nevertheless, a close 
inspection shows that the behavior of 1/t(u) is sublinear, i.e., 1/r oc u a where a < 1. 

Figure |3|b shows the comparable result for cold quasiparticles: we display the average 
scattering rate averaged over values of Ak, such that 0.5/£ < Ak < 2/£, for the same 
temperatures as in Fig. ||a. This average scattering rate is a good measure of the effective 
relaxation rate, l/r e //(u;), seen in optical measurements, since the behavior of quasiparticles 
at intermediate values of Ak, i.e., (Ak) max ~ 0.5 — 1, dominates the transport. Aside from 
numerical value, which is considerably reduced from that found at the hot spots (see Fig. 
|3]c), we see that the behavior is FL-like up to a frequency of order 15 meV, becoming linear 
with increasing u. There is no sign of sublinear behavior in 1/r. 

In Figure |]c, we have plotted the relaxation times for both hot and cold quasiparticles 
in YBa2Cu307 at T = 150K, in order to demonstrate the considerable variation in quasi- 
particle behavior as one moves around the Fermi surface: the frequency dependence of the 
hot and cold quasiparticle lifetimes are seen to be both qualitatively and quantitatively dif- 
ferent. Hot quasiparticles possess a mean free path of order of a few lattice spacings and 
are comparatively ill defined at higher u, while in the cold regions of the FS we find very 
well defined quasiparticles, with mean free paths of up to ten lattice spacings. Note that 
the mean free path is proportional to i>| and hence is rather sensitive to the choice of this 
parameter. Our choice oivf = 0.25eV is somewhat conservative, and the exact value of the 
mean free path is probably somewhat larger than that displayed in the Figure. 

We note that, due to the FL-like behavior at low frequency, the behavior of l/r(a?) 
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for cold quasiparticles at different temperatures becomes increasingly more distinct as T 
increases. This is in sharp contrast with the behavior observed in the singular (hot) region, 
especially when £ acquires the temperature dependence appropriate to the z = 1 scaling 
case. This is shown in Figure |]a, where we plot our result using spin fluctuation parameters 
appropriate for the underdoped YBa 2 Cu 4 8 material.0 The temperature dependence of 1/r 
is now completely suppressed at higher frequencies and eventually even inverted from the 
conventional FL behavior, especially at larger T. By contrast, for the cold quasiparticles, 
as shown in Fig. f|b, we obtain results which are very similar to those shown in Fig. ||b. In 
Sec. [IV] we will show that these results lead to the observed optical properties in HTSs.0 
In particular we show that 1/T e ff is linear in u regardless of the scaling law and hence, 
provided that the vertex corrections are not large, transport in underdoped materials above 
T* is very similar to that in optimally doped materials. 

IV. OPTICAL CONDUCTIVITY 

We now turn to the optical conductivity. We begin with the Kubo formula and the well- 
known expression for the optical conductivity in the absence of the vertex corrections:0 



<V(") = E« / duj'du"A k {u')A k {u") f{u ') f } U " ] (26) 



where N k is the number of points in the BZ and v k = Ve^, which in the case of a four- fold 
symmetric system, where x and y directions are equivalent, yields for the real part of a xx {uj): 

Rea xx (uj) = I duj'A k (uj')A k (uj' + uj) fiuJ,) ~ f{uJ ' + Uj) . (27) 



One can, of course, obtain Im a using the Kramers- Kronig relations. In Eqs. (|26| ) and 
Ak(cu) is the single particle spectral function, 

H ^ [^-e fe -i?e£ fe H] 2 + [/m£ fc M] 2 1 ] 

and £fc(u>) is the single particle self-energy. 
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In general, Eq. (p7|) is hard to solve analytically. However, depending on the proper- 
ties of the self-energy, various approximations can be used. For example, if the self en- 
ergy is only weakly momentum and frequency dependent, one can use the relaxation time 
approximat ion :0 

which is the result one obtains using the Boltzmann transport theory. It has been argued 
that qualitatively the optical properties of cuprates can be captured within this formalism, 
although, assuming a weak dependence of H(k,w) on k and u clearly represents a consider- 
able oversimplification. 

An interesting result can be obtained from Eq. (p9j). On performing the appropriate 
integrals, at low T, we find that 



2 c T , 

Rea xx = e — log 



1 + 



CT 



(30a) 



and 



Ima xx = e 2 ('^ + CT tan^ 1 - GT%) (30b) 



where C = 4vf/ag 2 . We see that, in contrast to ordinary FLs, in NAFLs the real part of a xx 
at high frequency depends only weakly on temperature. This weak temperature dependence 
of a xx (u>) at high frequency is a generic feature of any system with a low energy soft mode, 
be it a NAFL, as above, or a system near a charge density wave instability. Note that Eq. 
(j30|) reduces to the appropriate FL results in the u — > limit. 

A generalization of the relaxation time approximation, the memory function formalismjl! 
is appropriate for a self energy which, in the vicinity of the FS, does not depend on mo- 
mentum, as is the case when the dominant scattering mechanism is due to phonons, as in 
ordinary metallic superconductors; at sufficiently high frequency, it leads to a Drude-like 
conductivity: 

°xx{u) = —. f-r . u v (31) 
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with a frequency dependent effective scattering rate, and a frequency dependent mass m*, 
m, = m e (l + A (a;)). This expression is not directly applicable for the case of highly 
anisotropic scattering rates such as those studied here. However, we show below that in 
NAFLs, at sufficiently high frequencies, one can obtain an expression of this form, in which 
l/r e jj represents an average scattering rate as a function both of momentum and frequency. 

We proceed in the following way: We first recall that in the NAFL model the imaginary 
part of the self energy depends only weakly on the component of the wavevector k perpen- 
dicular to the FS; the dominant momentum dependence comes from the component which 
lies in a direction parallel to the FS (see Eq. (|19D and the discussion which follows). This is 
clearly an approximation, but not an unreasonable one, since only quasiparticles in a narrow 
region (of order T) around the FS are involved in the scattering (see Eq. fl27|)). In this case 
one can perform the change of variables (0) and integrate over e' explicitly. If the integral 
can be extended to infinity, the result can be expressed in the quite compact form: 

itaUoO = 4 / A id J^-^' + ^ lJ - j- (32) 

XXK ' 87T 2 JFs\v f \J UJ \u-E k (u+LU') + E* k (u')J V ' 

In the limit our > 1 we can expand this expression and compare our result to the Drude- 
formula in the same limit. We find: 



T e ff((j) JFS 2 i LO \T k (u + 0j') T k (u') 

In the limit u ^> T this yields 



(34) 



where the brackets indicate the appropriate average. We have already obtained the momen- 
tum averaged scattering rate, Eq. ([23|), which, at high frequencies, took the MFL form, Eq. 
( p5|) . These expressions show that the average (over momentum) of l/r k is a simple power 
law as a function of frequency, being quadratic and linear at low and high u respectively. As 
a result, averaging over frequency will not qualitatively change the frequency dependence 
of the momentum averaged scattering rate, < l/r k > k , since an average of a power-law 
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function is that same function. Therefore, aside from numerical factors (of order 2) and a 
somewhat different crossover scale E , the momentum averaged scattering rate < 1/V& >k 
is a reasonable approximation for l/r e //(a;) in the limit, wr e // 3> 1. 

Several important conclusions can be drawn from this result. First, only at high frequen- 
cies should the Drude formalism be used to analyze experiments. But, since at sufficiently 
high frequencies, in general, the quasiparticle picture breaks down, the results we obtain are 
likely only qualitative. One can show that Drude-like behavior is maintained as one goes to 
lower frequencies, but the corresponding scattering rate is no longer simply related to the 
quasiparticle self-energy, and will not be frequency independent. Moreover, quite impor- 
tantly, in the limit 1/ur <C 1 any feature seen in experiments, such as the weak temperature 
dependence of l/r e //(o;) at higher T and u>, or even the pseudogap seen in underdoped 
materials, must stem from the quasiparticle lifetimes analyzed in the previous section, since 
it is their average which is probed. 



V. COMPARISON WITH EXPERIMENT 

We now compare our results with experiment, in order to test the applicability of NAFL 
theory when one uses realistic parameters for both the effective interaction and band pa- 
rameters and carries out the relevant integrals numerically. Since the optical conductivity 
as a function of frequency frequently appears featureless, it is customary to use the mem- 
ory function formalismB to extract the effective transport scattering rates from experiments. 
Where this has been done, the analysis of the previous section allows us to compare directly 
our (averaged) scattering times with experiment. 

For systems for which the scattering rates have not been determined we compare our 
calculated values of a xx (u) directly to experiment. In so doing, we make use of our numerical 
results for the optical conductivity obtained by substituting the self energy (|5]) into Eq. 
(PT|), with and without including the inclusion of the vertex = cos 6 for I/tj. where 6 
is the angle between k and k! . Strictly speaking, this form of the vertex correction is only 
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applicable if the effective interaction is weakly frequency dependent, since only then is the 
current vertex proportional to the current, and hence additive to the self-energy. In NAFLs 
this is a somewhat questionable approximation, but as long as qvf <C 1, while V e ff is highly 
peaked near Q and the FS is large, the magnitude of the vertex correction turns out to be 
comparatively small, leading to corrections of order 10% for systems near optimal doping. 
However, since it is straightforward to include it numerically, we retain the function 7^/ in 
the subsequent numerical results. 

We begin with overdoped materials. While it is extremely important to compare the 
calculated optical conductivity with experiment only for the samples which are used to 
determine the spin fluctuation parameters, NMR and optical experimental results on the 
same sample are often not available. Thus for the first sample we consider, single layer Tl 
2201, the samples used by Puchkov et al,0 for optical experiments had a T c = 15K, while 
the best available NMR experiments on a comparable sample are those of Itoh et al,El on 
a Tl 2201 sample which had a T c = 23K. For overdoped systems, the various parameters 
(a, u s f, £) which determine the spin fluctuation spectrum do not vary appreciately with 
the doping level; hence our use of parameters chosen to fit the NMR experiments of Itoh 
et al should not produce an appreciable error. In choosing band structure parameters for 
our calculation, we assumed the band structure for the Tl system is close to that found in 
BSCCO family. Our calculated results, are compared with experiment in Fig. [5[ We present 
our results only for T > 100K, since for T <C 100K the stability of our numerical solution 
becomes questionable. Clearly the fit is quite good; both theory and experiment show that 
the system crosses over from the quadratic to linear in frequency behavior at a frequency 
7tT , which is of order uj s f$ i 2 (Ak)^ nax /2 = 40meV. Note that the crossover behavior changes 
very little with temperature, as predicted by Eq. (^) for the case of the z = 2 scaling present 
in all (sufficiently) overdoped materials. We further note that only two fitting parameters 
appear in our calculation: the coupling constant, g, and the zero frequency offset (which is 
very small, 1/t(u = 0) ~ 30 cm -1 and is likely generated by impurities). If the resistivity 
were measured for this material both fitting parameters would be fixed by experiment. 
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In Fig. H we compare our calculated optical conductivity for optimally doped YBa2Cus07 
with the experimental results of Ref. |20|. We have assumed, as in previous work, that 
uj s f = [6 + 0.06 T (K)]meV, ^ s /£ 2 = 60meV. Strictly speaking this scaling law, applicable 
only for temperatures above T cr ~ 140K, is only approximately correct at lower temper- 
atures. At high temperature our calculated results agree quite well with the experiment, 
both qualitatively and quantitatively. At low T qualitative agreement is found; the calcu- 
lated (7 n (w) at T = 100K crosses the calculated a xx {oj) at T = 300K at roughly the same 
frequency as that seen experimentally, although at high to the calculated values of u xx are 
considerably lower than those seen experimentally. We attribute this discrepancy, in part, 
to the possible presence of disorder in the experimental samples: for example, we found it 
necessary to use a higher value of the coupling constant g (g ~ 0.6eV) than that used in Ref. 
|9] (g ~ 0.5eV) to fit the resistivity. A small amount of, e.g., impurity scattering can bring 
about significant changes in transport and optical behavior at relatively low temperatures, 
while producing a negligible effect at higher T. 

An interesting feature of the YBa2C x O y family is the a-b plane anisotropy, which is seen 
in untwinned single crystals near optimal doping. [We have earlier provided a qualitative 
explanation for the anisotropy in the resistivity of these materials! and it is instructive to con- 
sider its optical counterpart.] In Fig. [7] we show the optical conductivity for a model system 
which displays strong a — b planar anisotropy, comparable to that observed in YB^C^Og. 
We follow Ref. [| and assume that the main influence of the CuO chains is to modify the 
electronic properties of the a-b plane and so model the observed anisotropy by introduc- 
ing anisotropic hopping matrix elements along the two different crystallographic directions, 
with t a = 0.554. This value of t a is somewhat low, but is not unreasonable for a qualitative 
analysis for pedagogical purposes. The input parameters (u) s f = [2 + 0.02 T (K)] meV and 
Lu s f£, 2 — 60 meV) are the same as in our earlier work.1^ As can be seen from the figure, the 
strong anisotropy we found for the d.c. resistivity has its analog in the case of <7 xx (u}), as 
might have been expected; the strongly conducting a direction displays a distinctly different 
uj dependence at higher frequencies. This behavior can be understood by recalling that 
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when the current runs along the a direction, the distribution of hot spots is such as to yield 
scattering rates which are mostly FL like, and hence one should expect to see a large con- 
ductivity even at large to. On the other hand, for current running along the b direction the 
dominant contribution to the conductivity comes from quasiparticles which possess strongly 
anomalous scattering rates; this leads to a much stronger frequency dependence of Rea^iuj), 
of the kind observed in many underdoped samples. 

We consider next the underdoped materials. Above T*, we expect the behavior of 1/Ve// 
to take the MFL form, Eq. (|25l). Indeed, Eq. (|25|) , obtained for a moderately large FS, 
does not depend on the spin fluctuation parameters, u s f and £, and hence independent of 
the scaling law relating u s f and £; 1/T e // is linear in both T and u> both above and below 
T cr , as long as T > T*. On the other hand, it is now widely believed that in heavily un- 
derdoped materials, below T*, the hot quasiparticles, which are close to the magnetic BZ 
boundary become effectively gapped (the pseudogap), with much of the quasiparticle weight 
transferred from low to high frequencies; the quasiparticle band dispersion thus acquires a 
form similar to that encountered with a preformed spin-density-wave (SDW) state. i Hence, 
at low frequencies, hot quasiparticles barely contribute to transport, a conclusion which is 
consistent with the superffuid density measurements at low temperatures, which suggest 
that only a fraction of the doped holes become part of the superconducting condensate. 00 
Although quasiparticles are not well defined in the hot region of the BZ, since their quasi- 
particle residua are rather small, they still contribute to transport up to frequencies of order 
T*. Above this scale their scattering is so strong that ImH is essentially temperature and 
frequency independent (see the result for Ak = in Fig. |2|). On the other hand, the cold 
quasiparticles are still well defined, albeit in a somewhat narrower region, centered around 
k = (n/2, tt/2), and symmetry related points in the Brillouin zone. However, since so much 
of the FS is gapped, the phase space available for the scattering of cold quasiparticles is 
considerably reduced. Thus for small frequency, lu < T*, one expects l/r e // oc u 2 , just 
as one finds 1/r ~ T 2 at zero frequency and temperatures well below T* in underdoped 
materials. In other words, there is a significant temperature variation of 1/T e // up to this 
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frequency. As the hot quasiparticles cease to contribute to transport above this energy scale, 
which is valid as long as the correlation length is large (that is, for £ > 2, and T < T cr ), 
the slope of l/r e // is reduced and it acquires its usual linear in u dependence. In terms 
of the NAFL model, a close study shows that the effective rate l/r e //(a;) increases up to 
frequencies of order T* ~ u> s f£ ) 2 (Ak)^ nax /2, where Ak max is much smaller (~ 0.1-0.5) than 
is found in optimally doped materials. This phenomenon is the optical analog of the re- 
sistivity which is a considerably stronger function of temperature below T* than above it. 
It is noteworthy that at sufficiently high temperature, above T cr , much of the quasiparticle 
weight is again transferred to the low frequencies and the hot quasiparticle again start to 
contribute to transport, but the slope of 1/V e // does not change, since it does not depend 
on the scaling law between the spin fluctuation parameters. 

We note that the fact that the contribution of the hot regions is strongly reflected in 
experiment in the underdoped materials suggests that large parts of the FS are in fact hot. 
Moreover, at higher frequency and low temperatures the slope of 1/T e // is the same as that 
found at high T in systems where there is no pseudogap present. If large parts of the FS 
exhibit highly anomalous behavior, then one should expect that at relatively low frequencies 
l/r e ff would begin to saturate. Indeed, we have found that the MFL form of l/r e ff, Eq. 
(p5j), is obtained only for a large FS, while for a smaller FS one would expect 1/T e // to display 
sublinear behavior as a function of frequency. This suggests that the linear in uj behavior 
of 1/T e // above reflects scattering processes involving momentum transfers far from Q, 
which supplement the magnetic scattering we have considered. If one is to maintain the linear 
behavior of 1/ r e f f for ttT and to larger than co> s /£ 2 , the dominant scattering mechanism must 
crossover from processes involving small to those involving larger characteristic frequencies. 
For example, in Ref. ^ we introduced FL scattering in order to avoid resistivity saturation 
at higher temperatures in 214 systems and obtained thereby quite reasonable agreement 
with experiment near optimal doping level. In this respect, optimally doped materials are 
merely those which have the smoothest transition from low energy to high energy scales, in 
excellent agreement with the experimental findings of Takagi et al.0 
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The additional scattering processes introduced to maintain linear in frequency 1/V e // 
above the pseudogap scale, do not qualitatively change the behavior of the effective re- 
laxation rate at low frequencies which is still governed by the strong magnetic scatter- 
ing. We present in Fig. [8] the result of a numerical calculation of l/r e //, where we have 
used Ak max = 0.4, taken u s f and £ to be the same as in Fig. [|a, and have incorpo- 
rated in our numerical work, based on Eqs. (p6|) and fl27|) , a FL-like scattering potential, 
Vfl = Xo/0- ~ ^^/r), where T = 400 meV,. Notice that l/r e // loses its temperature 
dependence at high frequency, i.e., all curves shown in the figure converge at a frequency 
~ u s f^ 2 (Ak max ) 2 /2. Below this frequency l/r e // displays considerable temperature depen- 
dence, even though we have not included the increase of u s f at low temperatures and a 
possible temperature variation in the effective size of the FS. Within the magnetic scenario, 
these scattering processes are closely related to the experimentally measured energy scales 
for q ps Q. Moreover, since the scale of uJ s f£, 2 is weakly doping dependent and always 
of the order of 1/tt times the Fermi energy the crossover energy should be almost mate- 
rial independent, in agreement with NMR experiments which show that for most materials 
£(T cr ) ps 2. 



VI. CONCLUSION 



In summary, we have applied both analytical and numerical techniques to the study 
of the transport properties of HTSs at optical frequencies within the NAFL model. We 
find that the peaking of the effective magnetic interaction leads to strong anisotropy of 
the quasiparticle properties in different regions of the BZ, and brings about a complex 
morphology for the optical conductivity and the effective scattering rate, l/r e //(o;), which 
can be extracted from it. We found that for a wide region of temperature and frequencies 
the NAFL model leads directly to the MFL form for l/r e //, Eq. fl2"5|), a scattering rate which 
agrees with experimental results found in nearly all optical measurements to date. Overall, 
in the appropriate limits, our results agree qualitatively and many cases quantitatively, with 
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the experimental findings. 

In this paper we have focused on the optimally doped and overdoped materials. In the 
underdoped case, to obtain results of comparable quality, one must parametrize in some 
detail changes in the quasiparticle behavior brought about by the quasiparticle pseudogap, 
which we have not yet done. However, to the extent that the FS undergoes a substantial 
transition as one varies doping from optimal to underdoped we have shown that the magnetic 
scenario provides results in qualitative agreement with experiment, even within this simple 
theory. 

Our use of the Born approximation, would seem to work somewhat better than one 
would expect it to in a system with strong correlations. Ideally, one should obtain H(k,u) 
self-consistently, before substituting it into Eq. fl27|). However, it may prove easier to cal- 
culate the current- current correlation function directly.il The results of calculations which 
include both the pseudogap and strong coupling effects will be presented in forthcoming 
publications.0'@ 
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FIGURES 

FIG. 1. A model FS in cuprates. For a quasiparticle interaction peaked at Q, the intersection 
of the magnetic Brilloiun zone (shown with the dashed line) with the FS specifies the singular 
points at which the interaction is maximally effective in low energy scattering processes. 

FIG. 2. The frequency dependence of the quasiparticle scattering rate at T = in (a) the 
case of a large FS. (b) the case of a small FS, for which Ak max = 0.3. In both panels the curves 
correspond to (from top to bottom) Ak = 0, 0.1, 0.25, 0.5, 1. The insets show the corresponding 
results at zero frequency as a function of temperature. 

FIG. 3. The frequency dependence of the quasiparticle inverse mean free path (in units of 
inverse lattice spacing), calculated for hvj = 0.25 eV, at (from bottom to top) T = 100, 200 
and 300K. The spin fluctuation parameters (uj s f = [6 + 0.06 T {K)\ meV, w s /£ 2 = 60meV) are 
appropriate for YBa2Cus07. Panel (a) shows the case with Ak = (hot region); panel (b) shows 
the momentum averaged inverse mean free path for cold quasiparticles (0.5/£ < Ak < 2); in 
panel (c) we compare the frequency dependences of this averaged cold quasiparticle mean free path 
(bottom) with that of hot quasiparticles (top) at T = 150K. Note the separability of the T and co 
contributions to the scattering rate in panels (a) and (b). 

FIG. 4. The quasiparticle inverse mean free path as a function of u at (from bottom to top, 
on the left hand side) T = 100, 200 and 300K in an underdoped material. The spin fluctuation 
parameters (u sf = [3.4 + 0.043 T (K)] meV, and £ = 50meV/u sf « 0.07 + 0.00086 T{K)) are 
appropriate for the pseudoscaling regime in a YP^C^Os sample. Panels (a) and (b) correspond 
to the hot and momentum averaged cold regions. Note the unusual temperature dependence of 
1/vfT at high frequency for the hot quasiparticles. 
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FIG. 5. A comparison of the calculated scattering rate as a function of co in Tl 2201 with the 
experimental results of Puchkov et al. 1 The curves correspond (from bottom to top) to T = 120K, 
200K and 300K respectively. The theoretical curves were obtained using spin fluctuation param- 
eters appropriate for this material, oj s f = [12 + 0.003 T (K)\ meV and uj s f£, 2 = 50meV. The FS 
size was estimated from the ARPES results on the closely related BSCCO compounds, for which 
Ak max srf 1.4. The coupling constant g 2 was adjusted to yield the correct spread of the curves at 
ijj = and a constant term was added, to take into account the likely presence of imperfections in 
the samples. 

FIG. 6. A comparison of the calculated optical conductivity of YBa2Cus07 with the experi- 
mental results of Tanner et al. 20 The top curve and the circles correspond to T = 100K and the 
bottom curve and the squares correspond to T = 300K. The parameters are those used earlier, 9 
uj sf £ 2 = 60 meV, u sf = [6 + 0.06 T (K)} meV, with t = 250meV and t' = -0.4t. Just as was the 
case for the resistivity, agreement with experiment is much better at higher temperatures. The 
value of the coupling constant g used here, g = 0.6eV, is somewhat larger than g = 0.53eV used in 
Ref. 9 to fit the resistivity, p xx . 

FIG. 7. Schematic behavior of a system with strong a-b plane anisotropy, such as YBa2Cu40s 
material. The lower (upper) curve corresponds to a aa (abb), with the solid (dashed) line referring 
to T = 100K (T = 300K). The input parameters (u sf = [2 + 0.02 T (K)] meV, uj sf ^ 2 = 60 meV, 
t a = 0.55*^) are the same as in Fig. 22 of Ref. SP. Clearly, both the magnitudes and the qualitative 
behavior of the two conductivities are different: abb is far more FL-like than is a aa . 
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FIG. 8. l/r e jf, obtained for a small FS, Ak rnax = 0.4, by numerically solving Eqs. (|i~7|) 
and (|26|), for the same input parameters as in Fig. 4, but adding a FL scattering contribution, 
Vfl = Xo/(l ~~ vkuj/T), with T = 400 meV, to the magnetic interaction, Eq. (4). The curves corre- 
spond (bottom to top) to T = 100, 150, 200, 250K. At high frequency the averaged quasiparticle 
lifetime is nearly T and oj independent, in accord with experiment on the underdoped materi- 
als. The stronger temperature dependence at low frequency stems from the strongly temperature 
dependent spin fluctuation spectrum found in all underdoped materials. 
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